Phase-locked magnetoconductance oscillations as a probe of Majorana edge states 
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We calculate the Andreev conductance of a superconducting ring interrupted by a flux-biased 
Josephson junction, searching for electrical signatures of circulating edge states. Two-dimensional 
pair potentials of spin-singlet <i-wave and spin-triplet p-wave symmetry support, respectively, (chi- 
ral) Dirac modes and (chiral or helical) Majorana modes. These produce /i/e-periodic magnetocon- 
ductance oscillations of amplitude ~ (e 2 /ft)A r ~ 1/ ' 2 , measured via an iV-mode point contact at the 
inner or outer perimeter of the grounded ring. For Dirac modes the oscillations in the two contacts 
are independent, while for an unpaired Majorana mode they are phase locked by a topological phase 
transition at the Josephson junction. 



I. INTRODUCTION 

Two-dimensional superconductors can support prop- 
agating edge states that are not localized by disorder 
for topological reasonsP^as a superconducting analogue 
of the metallic edge states in the quantum Hall effect 
or quantum spin-Hall effect Unlike the dispersionless 
"flat band" edge states of nodal superconductors P which 
leave a strong signature in the density of states, the prop- 
agating edge states have not yet been observed. They 
have been predicted in a variety of materials — includ- 
ing strontium ruthenateP heavily doped graphene,^ and 
topological insulators on a superconducting substrateP 

The symmetry-based classification of topological su- 
perconductors lists three types of propagating edge 
states.^ chiral Dirac modes, and chiral or helical Majo- 
rana modes (see Table |l| . A spin-singlet superconductor 
with d x 2_ y 2 + id xy orbital pairing supports edge states 
that propagate in one direction only (chiral) and are not 
sclfconjugate (Dirac fermions). For spin-triplet p x + ip y 
pairing the edge states are chiral and selfconjugate (Ma- 
jorana fermions). Counterpropagating (helical) Majo- 
rana modes are also stable against localizationpM31 un- 
like counterpropagating Dirac modes. 

The topological protection allows for correlations in 
the electrical current measured at distant points on the 
boundary connected by an edge statef^22l p or example, 
in an early study of this type, Law, Lee, and Ng con- 
sidered a superconducting disc deposited on the surface 
of a three-dimensional topological insulator!^ A chiral 
Majorana mode is confined to the perimeter of the disc, 
when the surface outside the superconductor is gapped 
by a ferromagnetP Two point contacts attached to the 



pair potential 


edge state 


symmetry class 


singlet, d x 2_ y 2 + id xy 


chiral Dirac 


C 


triplet, p x + ip y 


chiral Majorana 


D 


triplet, p x ± ip y 


helical Majorana 


Dili 



TABLE I: The three types of propagating edge states in a 
two-dimensional topological superconductor. 




FIG. 1: Left panel: Superconducting ring (S) containing a 
weak link (J), forming a flux-biased Josephson junction. A 
current can be injected into the grounded superconductor 
from a voltage-biased normal-metal (N) contact at the inner 
or outer perimeter. The Andreev conductance G n = I/V n is 
the ratio of the current-to-ground I and the applied voltage 
V n to contact n = 1,2 (with V = for the other contact). 
Right panel: Scattering processes at the outer perimeter, in- 
volving the coupling of a chiral edge state to 2JV incoming and 
2N outgoing electron-hole modes at the normal metal. This 
coupling introduces a dependence of the conductance on the 
phase difference <j) = (2e/ft)<& across the Josephson junction 
and on its topological quantum number a. 



perimeter measure a correlated current, mediated by the 
circulating edge state, and dependent on the number of 
magnetic vortices inside the disc. It is essential that the 
contacts share a boundary. If the disc would be replaced 
by a ring, with one contact at the inner and one at the 
outer perimeter, then there would be no correlations. 

Here we revisit this problem of edge-state mediated 
correlations in a superconducting ring, to show that the 
physics changes qualitatively if the ring contains a weak 
link forming a Josephson junction (see Fig. [lj. The 
weak link is a one-dimensional subsystem of the two- 
dimensional topological superconductor, with its own 
topological phase transition! 2 ^ Since magnetic flux $ can 
enter into the ring along the junction, there is no flux 
quantization and we can ask for the ^-dependence of the 
conductances G\ and Gi measured between the grounded 
superconductor and either the inner or the outer perime- 
ter. Dirac and Majorana modes both produce h/e- 



2 



periodic oscillations in the conductances, but only Majo- 
rana modes can correlate the oscillations in G\ and Gi- 
The mechanism by which the inner and outer perimeter 
communicate is a closing of the excitation gap when the 
Josephson junction undergoes a one-dimensional topolog- 
ical phase transition. The same conclusion was reached 
recently by Wieder, Zhang, and KaneP2 

The outline of the paper is as follows. In the next 
section we formulate the scattering problem of Andreev 
reflection from the perimeter of a superconducting ring, 
to obtain a general formula for the electrical conductance. 
The conductance G($, a) depends in general both on the 
enclosed flux $ and on the Z2 topological quantum num- 
ber a of the Josephson junction. A topological phase 
transition switches a between +1 and —1, resulting in a 



jump 5G of the conductance. In Sec. Ill we calculate the 
probability distribution P(8G) in an ensemble of disor- 
dered rings, using the method of random-matrix theory 
These are model-independent results, which take as input 
only the symmetry class of the topological superconduc- 
tor. We then turn, in Sec. |IV[ to specific model Hamilto- 
nians in each symmetry class. The numerical results for 
these models are discussed in Sec. [V] to arrive at a set of 
experimentally observable signatures of (1) the presence 
of circulating edge states and (2) their Majorana or Dirac 
fermionic nature. 



II. SCATTERING FORMULA FOR THE 
CONDUCTANCE 

We formulate the scattering problem for a supercon- 
ducting ring interrupted by a Josephson junction, enclos- 
ing a magnetic flux <E>. The ring is contacted at the inner 
or outer perimeter by a normal-metal electrode. Far away 
from any gap closings there is no transmission between 
the inner and outer perimeter, so we can treat these two 
scattering problems independently. We calculate the An- 
dreev conductance G between the normal metal (N) and 
the (grounded) superconductor (S). An /i/e-periodic flux 
dependence of the conductance serves as a signature of 
edge states circulating along the ring. 

As illustrated in Fig. [l] the NS interface is described 
by a scattering matrix S, with submatrices r-^ (reflection 
back into the normal metal), i c dg C ,N, ^N,cdge (transmis- 
sion from the normal metal into an edge state, and vice 
versa), and £ e dge (transmission along an edge state with- 
out entering the normal metal): 



S = 



r N ^N.edge 



The scattering matrix «Sj of the Josephson junction 
describes how the edge states return back to the NS in- 
terface after encircling the ring, 

a cdgc = Sj a cdgc- (2-4) 

Elimination of the edge state amplitudes gives the rela- 
tion a£j ut = i?a™, with the effective reflection matrix of 
the NS interface, 

R = ru + iN,cd g c(l — S , jicdgc)~ 1 <5'jicdgo,N- (2.5) 
The matrix R is unitary, RR) = 1, with electron and 



hole submatrices, 



R 



-*^ee ^ch 



(2.6) 



describing normal reflection (from electron to electron or 
from hole to hole) and Andreev reflection (from electron 
to hole or vice versa). The linear response conductance 
(in the zero-temperature limit) is given by 

G = Go Tr (1 - R ee Rl + R hc Rl) = 2G Tr R h X c , 

(2.7) 

with Go = e 2 /h the conductance quantum. It it conve- 
nient to rewrite this without reference to the submatrices, 



G/G = \ Tr (1 - Rt z R?t z ) , t z = 




(2.8) 



The Pauli matrix t z acts on the electron and hole degrees 
of freedom. 

The edge channels of a spin-triplet p-wave supercon- 
ductor are self-conjugate Majorana modes. It is then use- 
ful, following Refs. ll5lTOl to transform from the electron- 
hole basis to the Majorana basis, 



R^URU\ U 




(2.9) 



Electron-hole symmetry at the Fermi level requires that 
the scattering matrix elements are real in the Majorana 
basis, so i?t = R T — R" 1 . Because the Pauli matrix t z 
transforms into Ut z U^ = —t v , the conductance is given 
by 



G/G = \ Tr (1 - RT y R T T y ) , r y = 







(2.10) 



In what follows we will work in the electron-hole basis 



(2.1) 



(2.8) for spin-singlet d-wave pairing (when the modes are 
not self-conjugate) and in the Majorana basis (2.10) for 
spin-triplet p-wave pairing. 



Incoming and outgoing wave amplitudes at the NS inter- 
face are related by 



^N^N + ^N.edgeO, 



go "edge 



out j. in 1 i in 

"edge — <<edge,N"'N ' t edge"'edge- 



(2.2) 
(2.3) 



III. RANDOM-MATRIX THEORY 

The effect on the Andreev conductance of a topological 
phase transition at the Josephson junction can be ana- 
lyzed in a model-independent way by means of random- 
matrix theory. We will first do this for an unpaired chiral 
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Majorana mode and then for a pair of helical Majorana 
modes. The Josephson junction cannot undergo a topo- 
logical phase transition for chiral Dirac modes, and will 
generically not for an even number of chiral Majorana 
modes, so these two cases are not considered in this sec- 
tion. 



A. Chiral Majorana mode 

The conductance depends on the magnetic flux in a 
way which is restricted by electron-hole symmetry at the 
Fermi level. The restrictions are most severe for an un- 
paired chiral Majorana mode: The only phase shift al- 
lowed by electron-hole symmetry is ir (mod 2tt), so the 
conductance remains flux independent except when the 
enclosed flux is h/Ae (mod h/2e). Let us investigate this 
case in some detail. 

A single Majorana mode corresponds to scalars Sj 
and tedgej to a row vector i e dge,N, and to a column vec- 
tor tN.cdge- The contraction ^cdgcNTj/^dge n produces a 
scalar, which vanishes because r y is an antisymmetric 
matrix. This eliminates one term when we substitute 
Eq. (pTa) into Eq. (|2.10|). What remains is 



G/Gq = \ Tr (f - r^Tyr^Ty 



1 - S,t 



icdgc,N T "y''N' 7 J/^N,edge- (3-1) 

J ^cdgc 

Because Sj is an orthogonal matrix consisting of a sin- 
gle matrix element, it can only equal ±f . Including a tt 
phase shift from the winding around the ring, we define 



a = -Sj E {+1,-1} 



(3.2) 



as the topological quantum number of the Josephson 
junction. The effective reflection matrix R of the NS 
interface, constructed from Eq. (2.5), inherits this topo- 
logical quantum number, 



Dot R = a. 



(3.3) 



This follows from general considerations for a topological 
quantum number in symmetry class Dp^ but one can 
check it explicitly from Eq. (2.5). 



The conductance for a ring with an unpaired Majorana 
mode is flux independent — except at topological phase 
transitions, when a changes sign and the conductance 
jumps by an amount ±SG with 



SG = G(a = 1) - G(a = -1) 



Using Eq. (3.1 1 this can be written as 
2 



SG/Gq 



1 



t 2 

edge 



^edge , N 7y ?"*N T y , edge ■ 



(3.4) 



(3.5) 



For a disordered NS interface we may consider an en- 
semble of scattering matrices S, generated by varying the 
disorder realization. A simple choice is the circular real 
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FIG. 2: Probability distribution of the change SG in the 
conductance of the superconducting ring when the Joseph- 
son junction switches from topologically trivial to nontrivial. 
These are results from random-matrix theory in the circular 
real ensemble (CRE, symmetry class D), for an unpaired chi- 
ral Majorana mode circulating along the ring and N modes in 
the contact to the normal metal. The solid lines for N = 1,2 



are from Eq. ( 3.6 1, the histograms are obtained numerically by 
averaging the scattering matrix S uniformly over SO(2iV + 1), 
and the dashed curves are the large- TV Gaussian limit (3.91. 



ensemble (CRE) of class-D random-matrix theoryP^UM] 
for which S is uniformly distributed in SO(2iV + 1): The 
group of orthogonal (2iV + 1) x (27V + 1) matrices O with 
DetO = 1. The integer N counts the number of modes 
in the contact with the normal metal, including the spin 
degree of freedom. The factor of 2 in 2N + 1 accounts 
for the electron-hole degree of freedom and the +1 refers 
to the unpaired Majorana mode. 

The effective reflection matrix R, constructed from S 
via Eq. (2.5), inherits the uniform CRE distribution in 



Cv(2iV) — the set of 27V x 2N orthogonal matrices with 
determinant a. The uniformity of R £ O a (2N) is a con- 
sequence of the uniformity of S € SO(2iV + 1) because 
the transformation S ^ S ■ (U ® 1) with U E SO(2iV) 
maps R onto R ■ Uq without changing the determinant. 

If N = 1 or N = 2 the distribution of SG follows 
directly from the knowrP^ distribution P a (G) of the con- 
ductance in the CRE: In both these cases P_(G) = 8(G— 
2Go) is a delta-function distribution, so we may equate 
P{5G) = P+(G = SG + 2G ). Since P+(G) = 6(G) for 
N = 1 and uniform in [0, 4Gq] for N — 2, we arrive at 



P(6G) = S(5G + 2G ), for N = 1, 

P(SG) = 1/4G , -2G < G < 2G , for N = 2. 



(3.6) 



For N > 2 the knowledge of P a (G) is not sufficient to 
determine P(SG), but it can be determined directly from 
the uniform distribution of S in SO(2iV + 1). We have 
carried out this calculation numerically for small N, see 
Fig. § 

For large TV we can approximate the matrix elements of 
S by independent Gaussians, of zero mean and variance 
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1/2N. We define the unit vectors 



(l-^d g e) 1/2. 



C1 _ +2 \l/2 • 
I 1 ''edge J 



so that the conductance change (3.5 1 is given by 

2N 

SG/Gq = -2 u„w m (r N ) m „. 



(3.7) 



(3.8) 



In the large-iV Gaussian approximation, SG/Gq is 
the sum of Gaussians with zero mean and variance 
(2/N)(u n v m ) 2 , so its distribution is again a Gaussian 
with zero mean and variance 



2iV 



Var (<5G/G ) = ^ ^ ^ n v m f = ~, N > 1. (3.9) 



n,m— 1 



B. Helical Majorana modes 

Spin-triplet pairing with time-reversal symmetry can 
produce a pair of counterpropagating (helical) Majorana 
modes. This is symmetry class DHL In the Majorana 
basis the scattering matrix is orthogonal, as in class D, 
with the additional time-reversal symmetry condition^ 



S 



TyS T Ty. 



(3.10) 



This is equivalent to the requirement that the matrix 
product S = iTyS is both orthogonal (S^ = S T = S^ 1 ) 
and antisymmetric (S T — —S). The class-DIII random- 
matrix ensemble (T-CRE) is generated by drawing a ma- 
trix O from the CRE and constructing 



5' 



ITyS = O ■ ITy ■ O 



O e SO(4M + 2). (3.11) 



The channel number M = N/2 again only refers to the 
orbital degree of freedom, each mode having a twofold 
Kramers degeneracy. 

The Josephson junction breaks time-reversal symme- 
try, for <j) 5^ (mod 7r), so it may couple the two edge 
states and cause backscattering at the junction. In the 
simplest description (not made in the numerical calcula- 
tions of the next section) we neglect this coupling and 
set 



5j 



1 0' 
,0 1, 



=> Si = ir v S. 



-l(JT v 



(3.12) 




1 1 
5G / G{) 



FIG. 3: Probability distribution of the conductance change 
upon a topological phase transition in symmetry class Dili, 
involving a pair of helical Majorana modes circulating along 
the ring. The contact to the normal metal has M orbital 
modes, each with a twofold Kramers degeneracy. The results 
are from the time-reversally invariant circular real ensemble 
(T-CRE) of random-matrix theory. The solid line for M = 1 
is from Eq. (3.161, the histograms are obtained numerically, 



and the dashed curves are the large- N Gaussian limit (3.17l 



and also inherits the uniform distribution of the T-CRE: 



R — ir,.R = O ■ 



T , OeO CT (4M). 



(3.15) 



We seek the probability distribution P(5G) of the con- 
ductance change upon a topological phase transition in 
the T-CRE. For N = 2M = 2 the knowrPl probabil- 
ity distribution P„{G) of the conductance gives sufficient 
information, since P-(G) = 6(G - 4G ) =>> P(SG) = 
P + (G = SG + 4G ), resulting in 



P(SG/G ) 



1 



8y/l + SG/4G 
for N = 2. 



-4 < SG/Gq < 0, 

(3.16) 



The distribution P(SG) for M > 1 has been obtained 
by generating random matrices O uniformly in SO (4M + 
2) and then constructing S in the T-CRE via Eq. ^3.11 >. 
Results are shown in Fig. [3] For N ^> 1 we have again a 
Gaussian distribution with zero mean and variance 



Var(<5G/Go) = -, N>1, 



(3.17) 



four times larger than Eq. (3.9) because of the twofold 
Kramers degeneracy. 



The Pfaffian of 5j is the class-DIII topological 
invariant!^ 



a = -PfSj e {+1,-1}. 



(3.13) 



The effective reflection matrix R, constructed from S 



and Sj via Eq. (2.5), inherits this topological invariant, 
Pf(iT y R) = a, (3.14) 



IV. RESULTS FOR MODEL HAMILTONIANS 

The analytical considerations of the previous section 
rely only on the fundamental symmetries of the Hamil- 
tonian, without reference to a particular model. Here we 
present numerical results for model Hamiltonians in the 
various symmetry classes. 
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A. Chiral pair potentials 

We consider a two-dimensional superconductor in the 
x-y plane, with pair potential A dependent on the mo- 
mentum p — —ih(d x ,d y ). The Bogoliubov-De Gennes 
Hamiltonian, in the electron-hole basis, has the form 




(4.1) 



It contains the single-particle Hamiltonian 

H=(p- eA) 2 /2m cS ~ f i + U, (4.2) 

with A(r) the vector potential, m c g the effective mass, p 
the chemical potential, and U(r) an electrostatic disorder 
potential. (We conveniently set the electronic charge to 

The electron-hole symmetry relations are different for 
the different symmetry classes, 

A D = -A D => Hb = -T x n* D r x , (4.3a) 



A c = A c ^ H C = -r y n* C T y . 



(4.3b) 



As specific models we take in class D the spin-triplet 
chiral p-wave pairing 



Ad = Pp^A^),^ + ip y }, 



(4.4) 



with operators symmetrized by {a, b} — \{ab + ba). In 
class C we take the spin-singlet chiral d-wave pairing 



A c = ^2(p - eA) a M a0 (p + eA) , 

a,/3 



M(r) = P p 2 A(r) 



1 i 
i -1 



(4.5a) 
(4.5b) 



Both pair potentials properly produce a gauge invariant 
Bogoliubov-De Gennes Hamiltonianpl 

e- lXT 'H{A, A)e ixr - = h(a~ ^Vx, e~ 2iX A^ . (4.6) 

Since A D = A Q e 10 and A c = A e 2ie when A(r) = A , 
A = 0, and p = pf(cos 9, sin 6>), the magnitude of the gap 
is independent of the orientation. We expect that more 
general anisotropic models will give the same qualitative 
results — provided that the gap does not vanish in any 
direction. 

The ring has a weak link of length -Router — dinner j 
with i?i nne r and i? utcr the inner and outer radius of the 
ring. We assume that the ring is wide compared to the 
London penetration depth Al but narrow compared to 
the Josephson penetration depth Aj, 

-Router -Rinncr « Aj. (4.7) 

The first inequality ensures that the magnetic field is 
screened from the ring except at the weak link, along 



which a flux $ can enter. The second inequality prevents 
vortices to appear inside the weak link. The gauge in- 
variant phase difference across the weak link then has a 
uniform value 4> = (2e/fi)$. We will use a gauge with 
a real uniform order parameter A(r) = Ao and a delta- 
function vector potential 



§9(-y)6(x)x, 



(4.8) 



for a Josephson junction at x — (aligned along the 
negative y-axis). 



B. Helical pair potential 

We construct a model Hamiltonian with helical pairing 
from two time-reversed copies of the class-D chiral p-wave 
pairing, p x Hp y . Spin-orbit coupling of the Rashba form 
couples the spin-up p x + ip y sector with the spin-down 
Px — ipy sector, promoting the symmetry class from D to 
DHL 



The Bogoliubov-De Gennes Hamiltonian (4.1 ) contains 
the single-particle Hamiltonian 

Hdiu = [{p-eA) 2 /2m cS -p + U]<r + a so (p x a y -p y a x ), 

(4.9) 

where c x , & y , cr z are the Pauli matrices acting on the spin 
degree of freedom and ao is the corresponding unit ma- 
trix. The spin-orbit coupling strength is denoted by a so . 
The helical pair potential is given by 



Adiii = p F 1 {A(r),p 2; cr z + ip y a }. 



(4.10) 



The electron-hole symmetry requirement in class Dili 
is the same as in class D, cf. Eq. (4.3), 



^DIII — Aj^jjj => Hdiii — — T xH 



i-DIII'z- 



(4.11) 



For A = and real A the class-DIII Hamiltonian satisfies 
the time reversal symmetry 



"Rdiii — VvHt 



(4.12) 



C. Topological phase transition at the Josephson 
junction 

The phase transition in classes D and Dili is evi- 
denced by a closing of the excitation gap at the Joseph- 
son junction when <fi = it (m od 2t t). The gap closing 
and reopening is accompanie d 23 * 26 -* by a sign change of 
the topological quantum number a = Det R (in class D) 
or a — FfiTyR (in class Dili). In Fig. g we illustrate 
the gap closing for the chiral and helical p-wave pairings 
(|4~4]) and Oil). 



Away from <p = it the gap immediately opens in class 
D, while the gap closing persists for some range of <fi 
in class DHL This is a consequence of translational in- 
variance along the weak link, see App. |XJ Only the gap 
closing at <fi = ir is topologically protected. 
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chiral p-wave 




FIG. 4: Excitation spectrum of an infinitely long Joseph- 
son junction along the y-axis, for different values of the 
phase difference <f>, calculated numerically from the discretized 
Bogoliubov-De Gennes Hamiltonian ETJP^Panel a) is for the 
class-D chiral p-wave pair potential ( |4.4|| and panel b) for the 
class-DIII helical p-wave pair potential ( |4.f 1 1 . The closing of 
the excitation gap at <f> = n is topologically protected. 




chiral d-wavc 



symmetric 
asymmetric 




ky/k-p 

FIG. 5: Same as Fig. [4] for a class-C chiral d-wave pair 
potential. The two curves are for (f> = n, with and without 
±x symmetry of the Josephson junction along the j/-axis. The 
gap closing now has no topological protection, but requires a 
spatial symmetry. 



In class C there is no gap closing that is topologically 
protected. However, as explained in App. [XJ the com- 
bination of translational invariance along the y-axis and 
x i— > — x reflection symmetry allow for a gap closing at 
<f> = 7r (mod 2 tt). Wc show this in Fig. [5] for the chiral d- 
wave pairing (4.5). Disorder will break these symmetries 



D. Numerical results 

We have discretized the Bogoliubov-De Gennes Hamil- 
tonian on a square lattice (lattice constant a, hopping 
amplitude t, see App.[B]). The geometry is shown in Fig. 
[I] with a pair of normal-metal leads (width 30 a) attached 
to the inner and outer perimeter (radii 50 a and 100 a). 

The leads are modeled by setting A , A, a so , and U all 
equal to zero, at a chemical potential /iN = 0.5 1 for which 
there are M = 6 orbital modes. Each of these modes 
is spin-degenerate when coupled to the chiral d-wave or 
helical p-wave superconductor, and nondegenerate when 
coupled to the chiral p-wave superconductor. 

At the weak link the hopping matrix elements are re- 
duced such that the transmission probability per mode 
is < 0.1. Disorder in the superconductor is introduced 
by a random on-site potential U(x,y), uniformly dis- 
tributed in the interval ( — ?7disorder/2, ^disorder/2). We 
tOOk C/disordcr = 0.7 t. 

We solve the scattering problem numerically,^ to ob- 
tain the scattering matrix S for the electron (e) and hole 
(h) modes incident on the superconductor from the inner 
contact (1) and the outer contact (2). We then calculate 
the Andreev conductances G\ and G 2 from the outer or 
inner contact to ground, 



« i = J Tr (1 - S UAe Ste,u + S lhtl Ah,ie) > (4-13) 



G2 = ^ Tr (1 - S 2ei2e S 2 e, 2e + S 2h A,2e) • (4-14) 

This expression holds also at the gap closing, when there 
is a nonzero transmission probability between contacts 1 
and 2, under the assumption that contact 2 is grounded 
for the measurement of G\ and contact 1 is grounded for 
the measurement of G 2 . 

To probe the gap closing we also calculate the thermal 
conductance 



6h 

+ ^2h,lh^2h,lh + S 2e lh S 2 e,lh) ' 



(4-15) 



and remove the gap closing. 



which measures the thermal current between the inner 
and outer perimeter at temperatures To and To + ST. 

Results are shown in Fig. [6] for several disorder real- 
izations. The results for the electrical conductance (top 
row) will be discussed in the next section, in connection 
with experimental probes for Dirac or Majorana edge 
modes. 

The thermal conductance (bottom row) is not easily 
measured, but is included here because it illustrates in 
a striking way the significance of topological protection 
for a gap closing. The change in sign of the topological 
quantum number at the class-D phase transition results 
in a peak of the thermal conductance that is quantized^ 
in units of the thermal quantum 7r 2 fcgTo/6/i, see Fig. |6jl. 
The gap closing in class C has no topological protection, 
there is no sign change of a topological quantum number 
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FIG. 6: Electrical conductance (upper row) and thermal conductance (lower row) in the ring geometry of Fig. [I] as a function 
of the enclosed flux <E>, for the three pair potentials (|4.4|, (4.111, and (4.5 1. The red and blue curves are calculated numerically 
for two different realizations of the disorder potential The thermal conductance is measured between the inner and outer 
perimeter of the ring. The electrical conductance is measured either between the inner perimeter and ground (Gi , solid curves) 
or between the outer perimeter and ground (G2, dashed curves). 



and no quantized peak, see Fig. [6f. The class-DIII gap 
closing has topological protection (no backscattering) if 
it happens when the flux is a multiple of h/Ae, so time- 
reversal symmetry is preserved. Disorder leads to small 
displacements of the transition away from h/Ae, allowing 
for backscattering and resulting in a small deviation of 
the thermal conductance peak from the quantized value 
(see Fig. [6^). 



V. DISCUSSION 

The numerical results of Fig. |6^,,b,c illustrate how cir- 
culating edge states manifest themselves in the magne- 
toconductance of the ring. All three types of edge states 
introduce a flux dependence with a period of twice the 
superconducting flux quantum <£> = h/2e. The magne- 
toconductance oscillations are sample specific, depending 
on the disorder realization. The inner and outer perime- 
ter experience a different impurity potential and thus 
show a different magnetoconductance, but with the same 
h/e periodicity. A measurement of the fundamental fre- 
quency of the Fourier transformed magnetoconductance 
would be an unambiguous way to establish the presence 
of circulating edge states. 

The magnetoconductance contains additional informa- 
tion, it can identify unpaired (chiral or helical) Majorana 
modes. These produce jumps SG in the conductance 
when the flux is close to an odd multiple of $o/2, associ- 
ated with a topological phase transition at the Josephson 
junction. Both the sign and magnitude of 5G is disorder 
dependent and different at the inner and outer perime- 
ter, but the flux $ c at which the conductance jumps lines 
up. Notice that even when different disorder realizations 
cause a small shift in <I> C (compare red and blue curves in 
Fig. [6^) , the conductance at the inner and outer perime- 
ter jumps at precisely the same $ c (compare solid and 
dashed curves). This phase locking is a striking signa- 



ture of a topological phase transition at the Josephson 
junction. 

A measurement of the fundamental frequency com- 
ponent cos($/2$o + a ) °f the magnetoconductance at 
the inner and outer perimeter of the ring would there- 
fore show a random and uncorrelated phase a for Dirac 
modes, and a correlated phase peaked at (modulo it) 
for an unpaired Majorana mode. 

These magnetoconductance signatures of Dirac and 
Majorana edge states can be helpful in the ongoing search 
for topological superconductors. Recent attention has fo- 
cused on hybrid structures combining strong spin-orbit 
coupling with induced s-wave superconductivity, to pro- 
duce an effective chiral p-wave pairingPSHSSl A super- 
conducting ring deposited on a three-dimensional topo- 
logical insulator would need a magnetic barr ier along 
the perimeter to confine the edge statesP^ Alterna- 
tively, one might induce superconductivity in the two- 
dimensional electron gas of a semic onduc tor heterostruc- 
ture with strong spin-orbit coupling j 33 * 34 ! such as an InAs 
quantum well, and confine the edge states electrostati- 
cally by gate electrodes. 
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Appendix A: Gap closings due to spatial symmetries 

In symmetry classes D and Dili the closing of the ex- 
citation gap at the topological phase transition of the 



Josephson junction is topologically protected, meaning 
that disorder cannot open up the gap. However, in Fig. 
[5] we see a gap closing in symmetry class C, where the 
Josephson junction remains topologically trivial. Mo- 
rover, in Fig. [4}) we see that the gap closing in the helical 
p-wave junction persists over a range of (f>, rather than 
being limited to a single (j> as it is in class D. Both these 
features are due to spatial symmetries, as we now explain. 

Translational symmetry along the weak link (the y- 
axis) permits us to consider the parallel momentum 
k y = q as an external parameter. The Hamiltonian 'H(q) 
describes a zero-dimensional system which can undergo a 
topological phase transition as a function of the param- 
eter q in symmetry classes D and BDlP^l At this transi- 
tion a Z 2 topological quantum number changes sign, so 
to open up the gap requires, either, the breaking of a 
symmetry, or the merging of a pair of gap closings at a 
single value of q. 

So how do we arrive in class D or BDI when we start 
out from class Dili or class C? As pointed out in Ref. |35] 
in a different context, spatial symmetries can do this. 

Let us first show that T-L{q) is in class D for helical 
p-wave pairing. On the one hand, the electron-hole sym- 
metry relation (4.111 gives 



U(q) 



(Al) 



on the other hand, the helical p-wave pairing has the 
additional symmetry 



(t z <E) <J y )H(q)(T z <E) cr y ) = H{-q). 



(A2) 



Taking these two equations together we arrive at a sym- 
metry relation for %{q) at one single value of q, 



Sm(q) = -H(q)tt, n = (T y (g> a y )K, 



(A3) 



with JC the operator of complex conjugation. Because f2 
is an anti-unitary operator that squares to +1, this places 
in symmetry class D, with a topologically protected 
gap closing. Indeed, as we see in Fig. |3Jd, the pair of gap 
closings at <p — tt persist as <fi is increased, until the gaps 
merge at q = 0. 

Turning now to class C, we will show that H(q) is in 
class BDI for chiral d-wave pairing at <f> — n and if the 
electrostatic potential U(x) is ±x symmetric. Firstly, the 
class-C electron-hole symmetry relation reads 

H(q) = -T y H*(-q)Ty. (A4) 

Secondly, for (f> = tt and A = the Hamiltonian is real, 

H*(q)=H(q). (A5) 

Thirdly, the combination of U{x) = U(—x) and A(x) = 
— A(— x) at (f> — tt gives 



T z Vn{q) = H{-q)r z V, 



(A6) 



where V is the reflection operator (x y- > —x). Eqs. (A4) 
and ( A6 ) together give 



The anti-unitary operator fi' also squares to +1. The 
symmetries (A5) and (A7) place H(q) in class BDI, pro- 
vided that <fi = tt and the reflection symmetry is unbro- 
ken. This is consistent with what is seen in Fig. [5] The 
gap closing for chiral c?-wave pairing can be removed ei- 
ther by increasing 4> away from tt or by breaking the ±x 
symmetry of the weak link. 



Appendix B: Gauge invariant discretization of the 
Bogoliubov-De Gennes Hamiltonian 

The discretization of the Bogoliubov-De Gennes 
Hamiltonian (4.1 ) with a momentum dependent pair po- 



tential requires special care to ensure that the resulting 
tight-binding model is gauge invariant. We go through 
the steps in this Appendix. Following the established 
procedure of minimal coupling, we first discretize with- 
out a vector potential, then perform a gauge transforma- 
tion on the lattice, and finally replace the gradient of the 
gauge field by the vector potential. 

The discretization for A — is carried out by replacing 
the differential operators by symmetric finite differences, 

d x f(x)^~[f(x + a)-f(x-a)], (Bl) 
la 

to arrive at the tight-binding Hamiltonian 



t(n, to) 



t cc (n,m) t ch (n,m) 
t hc (n,m) t hh (n,m) 



(B2) 



The indices n, m label sites r n ,r m of a square lattice 
(lattice constant a). The diagonal elements n = m are 
the on-site energies and the off-diagonal elements n =/= m 
are the hopping amplitudes between sites n and m. 

The single-particle kinetic energy gives the electron- 
electron matrix elements, 

t cc (n,n) = 4t- n + U(r n ), t = h 2 /2m cS a 2 , (B3) 

—t for n, to nearest neighbours, 



otherwise, 



and the hole- hole matrix elements thh(n,m) = 
-t es (n,m). 

For the chiral rf-wave pair potential ( |4.5| ), still at A = 
0, we obtain the nearest-neighbor hopping amplitudes 



t ch (n±ax,n) = [A(r n ) + A(r„ ± ax)} , 



1 



t ch (n ± ay, n) = ^ [A(r„) + A(r„ ± ay)] 



(B4) 



the next-nearest-neighbor hopping amplitudes 



n'U(q) = -U{q)9!, Q' = t x VK. 



(A7) 



Aq 2 
±i 

v 



t eh _(n + ax ± ay, n) = [A(r„ + ax) + A(r„ ± ay)} 



t e h,(n — ax± ay, n) = — [A(r„ - ax) + A(r„ ± ay)} 



(B5) 



9 



and the on-site matrix elements 

t ch (n,n) = — =-[A(tv, + ax) + A(r„ - ax) 

2<T (B6) 

- A(r n + ay) - A(r n - ay)]. 

(We have defined q = fcpa.) These are all hopping am- 
plitudes from hole to electron. The hopping amplitudes 
from electron to hole are related by Hermiticity, 



This is the lattice analogue of Eq. (4.6). 



The effect on the electron-electron and hole-hole hop- 
ping amplitudes is the Peierls substitution,^ 



t hc (n,m) = t* ch (m,n). 



(B7) 



t QC (n, m) = t cc (n, m) exp ( i— j A ■ dl ] , 



(B9) 



We now introduce the vector potential A(r) 
-(h/e)Vx( r ) by means of the gauge transformation 



t(n, m) = e 



t(n, m)e 



A(r„) = e- 2 ^")A(r„). 



(B8) 



The line integral of the vector potential is taken along 
the lattice bond from site m to site n, and with this 
prescription the Peierls substitution can also be applied 
to vector potentials that do not derive from a gauge field. 



The transformed electron-hole matrix hopping amplitudes for chiral rf-wave pairing are given by 



i e h(n + ax ± ay, n) 
t e h(n — ax ± ay, n) 



t ch (n±ax,n) = ^ f 
t ch (n±ay,n) = ^ 
if /••>+<•**»« A dl 

■d n J n — ax 



e *f /" ±a * Adl A(r„) + A(r„ ± ax)e~^ ^ Adl 



(BIO) 



4<? 2 
±i r 
V 



A(r„ + ai)e-S ^ Adl + & ±af> A ' di A(r„ ± ay)e~^ ^ A ' dl 
A - dl A(r n - ax)e~^ /„"— A-dl + e *f ±nS ^<«A(r„ ± ay)e~^ ^ A dl 



t eh (n,n) 



2q 2 



Adl 



A(r„ + ai)e- 4 T In A dl + A( r „ - a^e"*^ /,! 
A(r„ + a»)e-^ ^ A dl - A(r„ - ay)e~^ ^ A-di 



(Bll) 
(B12) 



A similar calculation for the chiral p-wave pairing (4.4) gives the electron-hole hopping amplitudes 



t ch (n ± ax, n) = ^ [e l t /» ±B " A dl A(r„) + A(r„ ± a£) e - 4 t /« 



1 Adl 



t eh (n±ay,n) 



±1 



(B13) 



There are neither next-nearest-neighbor hoppings, nor 
on-site electron-hole matrix elements in this case. 



Notice that the discretized p-wave pair potential (B13 ) 



depends explicitly on the vecto r potential, while in the 
continuum representation (4.4) the vector potential en- 
ters only implicitly through A(r). All of this is required 
by gauge invariance. 

Finally, we give the corresponding expressions for the 
helical p-wave pairing (4.10). There is now a spin degree 
of freedom a an d the pair potential is diagonal taken from Eq. (B13) 



in that index. The electron-hole hopping amplitudes are 
given by 

f e t,ht(ra ± ax, n) = -t c ^(n ± ax, n) = t ch (n ± ax, n), 

i e t,hf(rc ± ay, n) = t ci . h i(n ± ay, n) = t ch (n ± ay, n), 

(B14) 



where matrix elements without spin indices should be 
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